Asymmetric Franck- Condon factors in suspended carbon nanotube quantum dots 
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Electronic states and vibrons in carbon nanotube quantum dots have in general different location 
and size. As a consequence, the conventional Anderson-Holstein model, coupling vibrons to the dot 
total charge only, may no longer be appropriated in general. Here we explicitly address the role of the 
spatial fluctuations of the electronic density, yielding space- dependent Franck- Condon factors. We 
discuss the consequent marked effects on transport which are compatible with recent measurements. 
This picture can be relevant for tunneling experiments in generic nano-electromechanical systems. 



PACS numbers: 73.23.-b; 85.85.+j; 78.32.-k 

Introduction — Advances in miniaturization paved the 
way to the fabrication of nanodevices in which molec- 
ular systems become active elements of circuits pQ. 
Tunneling of electrons through molecules leads to the 
excitation/de-excitation of quantized vibrational modes 
(vibrons) which have been experimentally observed in 
suspended carbon nanotubes (CNT) [2H5]. Their remark- 
able electronic and vibronic properties allowed for the 
observation of breathing [2] and stretching vibrons [3j 0] 
in recent transport experiments. 

In general vibrons couple both to the total dot charge and 
to the spatial fluctuations of the electron density. The 
latter received limited attention so far [6HH]. In most 
cases the Anderson-Holstein (AH) model [9j [10] has been 
employed, in which the vibron couples only to the to- 
tal charge. The AH model yields position-independent 
Franck- Condon (FC) factors pTj which strongly affect 
transport [T2HT4] . The predicted current suppression at 
low bias and the intensity of the vibrational sidebands 
have been confirmed in a recent experiment on suspended 
CNT quantum dots gj. 

In this paper we show that the effects of density fluctu- 
ations are crucial when the size and location of the dot 
and of the vibron do not coincide. They are indeed dra- 
matic when the vibron size L v is smaller than the dot size 
La: here, in sharp contrast with the AH model, position- 
dependent FC factors arise, possibly asymmetric on the 
dot tunneling barriers. This has profound consequences 
on the transport properties of the system. Only when 
L v > Ld, the total charge contribution is dominant and 
an effective AH model may be justified [T5] . 

Our predictions find an important confirmation in fur- 
ther measurements on the device considered in Ref. 0J 
A scanning electron microscope image, Fig. [T^i, shows 
the CNT connected to source (S) and drain (D) leads. 
A central suspended electrode (TG) acts as an electri- 
cally insulated top-gate, below which a quantum dot is 
formed (for more details see Ref. 4). Transport mea- 
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FIG. 1: (Color online) (a) Scanning electron microscope of 
the suspended CNT (arrows denote its position) connected to 
source (S) and drain (D) contacts. A top gate (TG) is also 
present. Scale bar: 200 nm. (b,c) Experimentally determined 
differential conductance G (units e 2 /h) as a function of the 
top gate voltage V g (units V) and bias V (units mV). (d) 
Schematic view of the coupled quantum dot-vibron system. 
The thick part represents the quantum dot and the wiggly 
line the vibron. 

surements have been performed in a pumped 4 He cryo- 
stat with a standard lock-in technique. The differential 
conductance G (Figs. [l]3,c) exhibits an almost perfect 
fourfold degeneracy in the Coulomb blockade diamonds 
and a rich structure of sidebands due to the excitation 
of stretching vibrons. The energy of electronic excited 
states measured on the Coulomb diamonds yields a dot 
size Ld ~ 240 nm, while the separation of vibrational 
subbands of about 0.8 meV yields L v « 60 nm <Ld |. 
A striking feature is the suppression of vibrational side- 
bands with negative slope as the gate voltage is varied. 
While in Fig.[]J, with V g in the regime analyzed in Ref.[4j 
sidebands with both slopes are present, in Fig. [TJd for 
V g < those with negative slope are completely absent. 
Here we show that this behaviour requires asymmetric 
FC factors at the tunneling barriers between the dot and 
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the leads. We stress that such a suppression cannot be 
obtained within the AH model, even assuming strongly 
asymmetric tunnel barriers. The case of Fig.[]J is on the 
other hand consistent with quasi-symmetric FC factors, 
in the spirit of the standard AH model. 
In addition, an alternating pattern of positive and nega- 
tive differential conductance (PDC/NDC) is observed in 
all the explored voltages ranges. This fact will be ex- 
plained in terms of a dynamical trapping of dot states 
induced by asymmetries in the tunnel barriers. 

CNT Dot-vibron model — As a model for our system, 
we consider a quantum dot confined between y\ = and 
?/2 = £d along the CNT and a vibron clamped at yo 
and yo + L v , with —L v <yo<Ld for a finite over- 
lap between the two systems (see Fig. We describe 
the CNT-quantum dot as a Luttinger liquid with two 
valleys n = ±1 and two spin channels a = ±1 [16] 
employing standard bosonization techniques with open 
boundaries [T71 [18] (i.e. the electronic field satisfies 
^77,0- (0) = V^cr (^d) = 0). The bosonization picture is 
not essential in our analysis, but it simplifies consider- 
ably the formal treatment of the electron-vibron cou- 
pling. The dot Hamiltonian is composed of three terms 
H d = + + (h = 1, /i € {c+, c-, s+, s-}) 

= f (N c+ - N s f + ^ (Nl + N! + + Nl) , 

H « } = 5 EE . 

fi q=l 

The term describes the energy of N c + electrons in 
the dot for a given configuration with electrons with 
spin o~ in branch 77. Here, total (+) and relative (-) charge 
(c) and spin (s) modes have been introduced [16], with 
N c+ = £„ CT iV, N c . = £„ CT T]N V a, N s+ = E„ ff ^V 
and N s - = ^2 rj(T 'no'N r}(T . In addition, N g oc V g is the 
charge induced by the top-gate voltage V g1 E c is the 
charging energy and the Fermi velocity [19 . Collec- 
tive charge and spin excitations are described as bosonic 
modes in . The generalized position and momentum 
of mode \i are respectively x^ jq and p^^ with frequency 
^fi,q = KV^q/La and group velocity [19 . Finally, ' 
models a shift between the energy of the two valleys [20] . 
The lowest stretching vibron is described by the har- 
monic Hamiltonian H v — p\)2M + MujqXq/2, where 
M is the vibron mass, Uq = ttv s /L v its frequency and 
v s the stretching mode velocity [19]. Here, xo is the 
amplitude of the lowest vibron, with distortion field 
u(y) = V2xo sin [ir(y — yo)/L v ] along the CNT, and po 
is the conjugate momentum. In a CNT, v s < and the 
experimental estimates yield ujq < uj^^i [3, 4 . 



Electrons and vibrations are microscopically coupled via 

min[L d , 2 y +L v ] 

H d - V =c J dy[p^ + \y)+p^ + \-y)]d y u(y), (1) 

max[0,y ] 

where c is the deformation potential coupling con- 
stant pjJl (2TJ [22] and we have introduced the to- 
tal electron density of right movers p^\y) = 

Ttw^R^aiyWRrtAy) with ^RmAv) their Fermi op- 
erator [23l [24]. Notice that, while vibrations cou- 
ple to the c+ mode only, all four collective electronic 
modes are important for transport. In bosonized form, 
one has p£ + \y) = (N c+ /2L d ) + (l/2ir)d y c+ (y) with 
My) = E q> o e-^ 2Ld [e— *y/ L * (x c ^ q - 

ioJ~+ q Pc+, q ) + h.c] and £ the short wavelength cutoff. 
This expression of the density neglects the fast oscillat- 
ing terms due to mixed right and left-moving fermion 
fields and is reliable in the large total charge 7V C+ regime 
with 7V C+ ^> Lji\L N . This condition is experimentally 
satisfied in all the ranges of parameters analyzed in this 
paper. The coupling Eq. (HI can thus be decomposed into 

^d-t = c o^o^V c + and nf^ v = ^oV^E^li c q x c+, q , due 
to zero modes and plasmons, respectively. The lengthy 
but straightforward expressions of Co and c q will be de- 
ferred to a future publication [25]. We point out that 
Eq. accounts for the coupling between vibron and 
density fluctuations i^-v' ne gl ec ted in the AH model. 
The total Hamiltonian H = H d + H v + H d _ v is thus 
quadratic in the generalized coordinates and is diagonal- 
ized p~3j[26] (details will be given elsewhere [25 ) into 

#o = f (^c + - N g f + ^ [Nl +iV s 2 + +JV?_] + tff 

+ ^^4^+ ^Uw b l,v b W' ( 2 ) 

The sectors with /i/c+ are clearly unaffected by Eq. ([I]) . 
On the contrary, in the c+ sector new modes, created by 
al with energies Vt v emerge. For v > 1 they represent 
new collective electron modes (dressed plasmons), while 
for v = a vibronic excitation dressed by plasmons is 
obtained. The latter is the low-energy vibrational mode 
observed in the experiments. The energies Vt v satisfy 
^ = "1 + ET=i c V (^-«c 2 J. with % < co and 
fl v > ^ c +,v for ^ > 1 always. Note that we have reab- 
sorbed a polaron shift into E c [9]. 

Local FC factors — We can now study how the 
bosonized Fermi field ^r^Av) 02] is affected by the 
transformation above. As we study tunneling at energies 
smaller than the collective charge and spin excitations of 
the dot, we restrict the Hilbert space to the v = mode 
of the sector c+ only. Due to Eq. ([!]), the vibron oper- 
ators a v appear in the electronic field, whose truncated 
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form after the diagonalization reads [27] 
1>R,vAv) ~ ^^e- [AN+A - (j/)1[o o- aol e a + (y)[a o +001 , (3) 

V27T^ 

where Xr?,cr decreases N Vi(7 by one, A at = cq/ \/2MVLq and 



A±(y) = 



9=1 



c q F±(y) 



(4) 



With K = l + E^l^/("g 



c+ 



? ) 2 and F ± (y) = 



sin (nqy/Ld + 7r/4±7r/4). Note that both A at and A±(y) 
depend on the CNT and dot parameters and position 
only via yo, the length ratio S = L v /Ld, the ve- 
locities ratio v c +/v s , and the dimensionless coupling 
A m = c/(v 8 y/Mu^) [19 . The local FC factors [9j [13] 

X^) = 2^|(^-1,^r,^(^)|^,0| 2 describ- 
ing tunneling of an electron off the dot while changing 
the vibron number from I to V (I < V) have the form 



x w (y) 



-A 2 (y) 



IHv)} 



2(V- l) ^-l {xHy))] 2 (5) 



with X 2 (y) = [Xn + A_(?/)] 2 + A+(y) a position- dependent 
effective coupling and ^(x) the generalized Laguerre 
polynomials. This is the main result of our paper. 
The position dependence is entirely due to the cou- 
pling between the vibron and the density fluctuations, 
neglected by the AH model which instead predicts 
position-independent FC factors, with constant interac- 
tion strength Xn- When max[A±(?/)] ^> Xn the position 
dependence cannot be neglected, and the AH model be- 
comes questionable. This occurs for S = L v /L^ < 1 
(which is the case of our experiment) and a vibron lo- 
cated inside the dot: in this case indeed Aat = 0. Fig. [2k 
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FIG. 2: (a) Plot of X(y) for 5 = L v /L d = 0.1 and different 
positions of the vibron center y c = yo + L v /2: (thick solid) 
y c — — L v /4; (thick dashed) y c = 0; (thin solid) y c = L v /2; 
(thin dashed) y c = La/4; (thin dotted) y c = La/2, (b) Plot 
of X(y a ) vs. S (a = 1, 2) for y = and a = 1 (solid) ; a = 2 
(dashed). Notice the strong asymmetry for (5 < 1 and the 
symmetric A's for S = 1. Here, v c +/v s = 32 and A m = 3 (for 
a CNT waist ~ 1 nm) [15] . 

shows A (y) for S < 1 and different locations of the vi- 
bron. When the latter sits inside the dot (thin lines, 
for < yo < — L v ), X(y) is sizeable only in the vi- 
bron region. For vibrons partially outside the dot (thick 
lines), A at ^ and the position dependence of X(y) is 



weaker. For S > 1 (not shown), A at ^> X±(y) which im- 
plies X(y) ~ A at, and the spatially-independent FC factors 
of the AH model are obtained [T5] . 

Of particular relevance for transport is the value of the 
coupling at the position of the tunneling barriers, X(yi) 
and A (1/2). For S < 1 and a vibron located asymmetri- 
cally with respect to the dot center, they become very 
asymmetric (see the thin solid line of Fig. 2i), yield- 
ing strongly asymmetric FC factors. In Fig. 2d, X(yi^) 
are shown as a function of S < 1 for a vibron located 
near the left barrier. The couplings are strongly barrier- 
dependent and vibrational excitations are strongly sup- 
pressed for tunneling on the right. In the symmetric case 
5 = 1, dot and vibron occupy the same region of space 
and A (2/1) = A (2/2) ©• Notice however that A at = 0. 
The maximum value of the coupling for S < 1 is cru- 
cially sensitive to the ratio v s /v c + and the value of 
A m . The coupling of the dot to the breathing mode re- 
duces v c+ [3 [28], increasing v s /v c + and allowing to reach 
X(y 1 ) > 1 with X(y 2 ) <C 1. In parallel, recent measure- 
ments in graphene [29] report a large deformation poten- 
tial, which further increases A m . 

Transport properties — In order to address the elec- 
tronic transport we introduce the tunneling Hamiltonian 
coupling the dot to the leads (represented by the CNT 
portions outside the dot) 

H t = 5Z^^^R^,a(^)^R,r7,a(^) +h.C. , 

a— 1,2 77, cr 

where t aiV are tunneling amplitudes and ^R^^iVa) * s 
the right movers field for lead a . In sequential tunneling, 
transition rates are evaluated between eigenstates of Ho - 
Eq. ([2]). For tunneling into the state rj of the dot through 
the barrier a one has [TJl [30] 



p(in) 



'1*2 



+ 11 



; X lv {y a )f[AE+{-ir +l eV/2] 



where To = 27i\D|t2,+i| 2 /£ 2 and V is the leads density of 
states, while f(E) is the Fermi function with AE the en- 
ergy difference between final and initial dot states. Sim- 
ilar expressions hold for tunnel-out processes. 
The experiment allows to estimate the relevant param- 
eters: E c « 4.5 meV (via Coulomb diamonds), the av- 
erage To ~ 1 /xeV (via current traces), ^0 ~ 800/ieV 
(average vibron sideband separation) and k^T « 90/ieV 
(for T « IK). Since k B T ^> T the sequential regime 
is justified, Q ^ ^bT allows to resolve vibronic exci- 
tations while Q ^> T justifies a rate equation [14] [31] 
neglecting vibronic coherences [32]. The extremely rich 
scenario obtained for different asymmetries of left/right 
tunnel barriers A = \ti y7] \ 2 /\t2 y7] \ 2 and of the coupling be- 
tween leads and the two valleys 7= |£a,-i| 2 /l^*,+i| 2 will 
be discussed in detail elsewhere [25] . 
Here we focus on the relevant case to address the exper- 
imental data in Fig. [I] For V g < in Fig. [lj}, we found 
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FIG. 3: (Color online) Plots of the numerical differential 
conductance G (units e 2 /h) as a function of N g — N g — 
37TV F /2E c L d and V (units meV). (a) Density plot for A = 
1/20, 7 = 20, X 2 (yi) = 2.4, X 2 (y 2 ) = 0.1; (b) same as in 
(a) but for X 2 (yi) = X 2 {y 2 ) — 2.4; (c) same as in (a) but 
for A = 1/5, A 2 (j/i) = 1.8 and X 2 (y 2 ) = 0.6. In all panels, 
= 0.8 meV, k B T = 0.1 fi , E c = 4.5 meV, Ae = 0.48 meV 
and To = 0.8 fieV. For simplicity, only one resonance is 
shown. 

that the only possible parameter range compatible with 
experimental data is: X(y2) <C A(yi) , A < 1, 7 > 1 and 
Ae > kBT. The asymmetry of the FC factors is respon- 
sible for the strong suppression of negative-sloped side- 
bands, as clearly shown in Fig.[3^i. We want to stress that 
the absence of traces with negative slope is not achiev- 
able in the standard AH model with symmetric FC fac- 
tors, even in the presence of a quite strong asymmetry of 
the tunneling barriers, as shown in Fig. [3J3- This proves 
the need to go beyond the AH model. 
The alternating PDC/NDC traces can be addressed in 
our model by the three remaining constraints (on A, 7 
and Ae). The NDC is due to the creation of a bottleneck 
in transport: when 7 > 1 tunneling into states rj = +1 
is strongly suppressed leading to a dynamical trapping 
and NDC [33j [34], while the states 77 = —1 provides a 
fast pathway with ensuing PDC. A shift of the two val- 
leys Ae > is necessary in order to resolve the two 
channels. Finally, the asymmetry A < 1 allows to obtain 
the PDC/NDC pattern in all voltage regimes V ^ 0. 
Analyzing the experimental data, we observe that for 
V g > the suppression of conductance traces becomes 
less severe (see Fig. lib), suggesting more symmetrical 
FC factors as in Fig. [3]3, in line with the standard AH 
model. In this case NDC traces with both positive and 
negative slopes occur for V > 0, pointing at an asym- 
metry A weaker than in Fig. |3^i. The ultimate reason 
for the relative shift of electronic vs vibronic wavefunc- 
tions at different V g lies in the unknown details of the 
electronic and mechanical confinements. Our predictions 
could stimulate further developments of experimental se- 
tups with full control over these delicate aspects. 

Conclusions — Recent experimental data show the 
need of a theory beyond the usual Anderson-Holstein 
model of quantum transport in nano-electromechanical 
systems. Here we investigate this new issue by consider- 



ing the combined role of the electronic charge and density 
fluctuations in the coupling to mechanical deformations 
for suspended CNT quantum dots. When vibrons are 
asymmetrically embedded into a larger dot, position de- 
pendent Franck- Condon factors arise. The consequent 
marked effects in the transport characteristics allow to 
address experimental features which could not be cap- 
tured by the standard AH model. Our analysis can be 
easily extended to consider e.g. planar metallic contacts 
or the tunneling from a localized tip. For small vibrons 
embedded in larger dots a spatially-resolved injection of 
electrons would show a tunneling suppression sensitive 
to the vibron location, making our theory relevant for 
spatially-resolved scanning tunnelling microscope mea- 
surements as well. Similar effects could be expected also 
in systems of higher dimensionality, such as e.g. quan- 
tum dots embedded into suspended graphene sheets. 
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